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Research  Goals: 

The  long-range  goal  of  our  research  program  is  to  understand, 
using  laboratory  experiments,  numerical  modeling  and 
theoretical  analyses,  the  small-scale  mixing  processes  that  occur 
in  oceanic  stratified  layers. 


Objectives : 

The  objectives  of  research  carried  out  during  the  past  year  were 
to:  (i)  understand  entrainment  processes  across  density  interfaces 
subjected  to  differential  turbulent  forcing;  and  (ii)  develop 
theoretical  and  laboratory-based  descriptions  for  internal  waves 
and  mixing  processes  resulting  from  the  interaction  between 
turbulence  and  stratification.  The  first  part  of  the  study  is 
applicable  to  the  evolution  of  density  step  structures  observed  in 
the  atmosphere  and  oceans,  which  consist  of  turbulent  layers 
separated  by  density  interfaces  (Padman  &  Dillon  1989,  J.  Phys. 
Oceanogr.,  18,  1458-1462  ).  The  second  investigation  is  aimed 
toward  parameterization  of  processes  occurring  at  the  oceanic 
convective  layers,  for  example,  the  night-time  boundary  layer 
where  cooling-induced  convection  plays  a  significant  role  in  the 
generation  of  turbulence  (McPhaden  &  Peters  1992,  J.  Phys. 
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Oceanogr.,  22,  1317-1329)  . 

Approach: 

The  study  on  differential  mixing  was  mainly  laboratory 
experimental,  supplemented  by  the  development  a 
phenomenological  model  based  on  prior  studies  related  to 
entrainment  from  one  side  of  the  interface.  The  experiments 
were  carried  out  by  using  a  two-layer  stratified  fluid  system  of 
which  the  homogeneous  layers  were  agitated  by  oscillating . grids 
that  generate  unequal  turbulent  intensities;  the  interfacial 
migration  rate  and  the  properties  of  turbulent  layers  were 
measured.  In  the  theoretical  study  on  mixing  and  internal  waves 
the  techniques  of  rapid-distortion  theory  were  applied  to 
calculate  turbulent  parameters  (e.g.,  r.m.s.  velocity  and  other 
integral  properties)  in  and  around  a  density  interface  forced  by 
contiguous  turbulent  layers.  Laboratory  experiments  were  also 
carried  out  by  subjecting  density  interfaces  to  grid-induced 
turbulence  and  the  results  were  compared  with  the  model 


predictions . 
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Tasks  Completed: 

In  the  first  part  of  the  study,  the  dependencies  of  the  interfacial 
migration  rate  and  buoyancy  variations  in  the  layers  (and  hence 
the  buoyancy  fluxes)  were  determined  as  a  function  of  the 
governing  system  parameters.  A  simple  model  was  developed  to 
predict  these  parameters.  In  the  second  part  of  the  study,  a 
theoretical  model  was  developed  for  the  evolution  of  linear  and 
non-linear  waves  on  a  thin  thermocline  (density  interface) 
subjected  to  shear-free  turbulence.  The  predictions  of  this  model 
included  various  turbulence  quantities  and  the  interfacial  mixing 
rate.  State-of-the-art  experimental  techniques  such -as  laser-- 


Doppler  velocimetry  and  laser-induced  fluorescence  were  used  to 
measure  the  characteristics  of  the  flow  fields  associated  with  thick 
and  thin  thermoclines . 


Scientific  Results: 

The  experimental  results  on  interfacial  migration _ rates  in  the 
presence  of  differential  turbulence  forcing  were  in  general 
agreement  with  the  phenomenological  model  developed  as  a  part 
of  the  study  (Figure  1) .  This  model  is  based  on  the  notion  that 
the  net  interfacial  migration  velocity  is  the  same  as  the  resultant 
of  entrainment  velocities  that  occur  when  the  opposite  layer  is 
non  turbulent.  Some  discrepancies  were  noted,  however,  for 
which  plausible  explanations  were  provided.  The  interfacial 
migration  occurs  in  a  direction  that  opposes  the  gradient  of 
inertia  forces  of  the  integral-scale  turbulent  eddies  at  the 
interface;  it  continues  until  the  inertia  forces  on  both  sides  have 
the  same  magnitude,  whence  the  interface  becomes  stationary. 

Yet,  the  buoyancy  transport  though  the  interface  continues  and 
finally  the  system  runs  down  to  a  homogeneous  state.  The 
results  were  recast,  based  on  the  properties  of  mixed-layer 
turbulence  and  density  interface,  into  a  form  that  is  convenient 
for  direct  oceanic  applications. 


The  model  developed  in  the  second  part  of  the  study  was  used  to 
calculate  the  r.m.s.  velocity  fluctuations  in  and  around  a  thin 
density  interface,  the  r.m.s.  amplitude  of  interfacial  waves  and 
the  diapycnal-mixing  rate.  The  agreement  between  the 
predictions  and  previous  laboratory  results  was  good.  The 
present  laboratory  experiments  with  a  thick  density  interface 
(sandwiched  between  a  turbulent  layer  and  a  weakly  stratified  or 
an  unstratified  layer)  revealed  some  interesting  phenomena;  a 
fourth,  strongly  stratified,  layer  was  found  to  develop  in  the 
region  between  the  thermocline  and  the  turbulent  layer,  and  most 
of  the  wave-turbulence  interactions  were  confined  to  this  layer . 
It  had  a  velocity  maximum,  indicating  the  presence  of  many 
resonant  modes  of  waves  (Figure  2) .  Detailed  measurements  of 
the  vertical  velocity  structure,  internal-wave  parameters  and 
mixing  characteristics  pertinent  to  the  density  interface  were 
made  and,  whenever  possible,  the  results  were  compared  with 
available  theoretical  predictions  and  geophysical  observations 


Accomplishments: 

A  significant  accomplishment  is  the  demonstration  of  the 
sensitivity  of  diapycnal  mixing  to  the  interfacial  Richardson 
number,  in  addition  to  the  commonly  used  bulk  Richardson 
number.  The  number  of  possible  resonant  internal-wave  modes 
varies  with  the  magnitude  of  the  interfacial  Richardson  number 
and  the  energetics  of  mixing  are  dependent  on  the . structure  of 
internal  waves.  Also,  we  were  able  to  obtain  a  simple  bulk 
parameterization  for  entrainment  in  the  presence  of  differential 
turbulent  forcing,  which  can  be  applied  in  the  studies  of 
oceanic  thermohaline  staircase  structures,  for  example,  those 
described  by  Kelley  (1987,  J.  Phys .  Oceanogr . ,  17,  1633-1639). 
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ocean  and  its  subsequent  mixing  by  small-scale  turbulence. 
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Figure  1:  A  comparison  of  ihc  model  prediction  (dark  line)  with  the  experimental  results  (symbols),  for  the  case  of 
differential  turbulent  forcing  on  a  density  interface.  The  variation  of  the  non-dimensional  interfacial  position 
(Dl/DO)  with  non-dimensional  time  (f)  is  shown.  Here  the  ratio  of  eddy  diffusivities  of  the  upper  and  lower  layers 
is  0.5;  Dl  represents  the  depth  of  the  upper  layer  and  DO  is  its  initial  value.  RiO  is  an  initial  bulk  Richardson 
number.  The  constants  Al  and  n  are  based  on  the  case  where  turbulence  is  present  only  on  one  side  of  the  interface. 


Figure  2:  The  variation  of  the  normalized  r.m.s.  vertical  velocity  within  a  (thick)  density  interface,  for  the  case 
where  the  intcrfacial  Richardson  number  (based  on  the  intcrfacial  thickness)  Rij  >  rc~.  A  velocity  maximum  can  be 
seen  within  the  interface.  In  the  case  cf  a  thin  interface  with  Rij  <  7i~,  only  the  first  mode  of  internal  waves  is 
present  and  the  velocity  profile  is  flat.  This  shows  that,  in  calculating  diapycna!  fluxes,  Rij  should  be  considered 
along  with  the  usual  bulk  Richardson  number  Ri.  p  is  indicative  of  the  vertical  position  within  the  interface. 


